Entanglement and Subsystem Particle Numbers in Free Fermion Systems 



O 
(N 



Y. F. Zhangi, Huichao Li^, L. ShengiEl R- Shen\ and D. Y. Xing^ 

^National Laboratory of Solid State Micro structures and Department of Physics, Nanjing University, Nanjing 210093, China 

(Dated: May 18, 2012) 

We study the relationship between bipartite entanglement and subsystem particle number in a 
half-filled free fermion system. It is proposed that the spin-projected particle numbers can distin- 
guish the quantum spin Hall state from other states, and be linked to the topological invariant of 
the system. It is also shown that the subsystem particle number fluctuation displays behavior very 
similar to the entanglement entropy. It provides a lower-bound estimation for the entanglement 
entropy, which can be measured experimentally. 
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I. INTRODUCTION 

Topological phases of matter are usually distinguished 
by using some global topological properties, such as 
topological invariants and topologically protected gapless 
edge modes, rather than certain local order parameters. 
The integer quantum Hall effect^, fractional quantum 
Hall effect^, and band Chern insulators'^ can be charac- 
terized by Chern numbers or Berry phases^. The quan- 
tum spin Hall (QSH) effect-i^ and the three-dimensional 
topological insulators'" are characterized by the Z2 in- 
variant^ or spin Chern number— lii. Recently, quantum 
entanglement^'^, 'which reveals the phase information of 
the quantum-mechanical ground-state 'wavefunction, has 
been used as a tool to characterize the topological phases. 
As sho'wn by Levin and We\^ and also by Kitaev and 
PreskiU^^, the existence of topological entanglement en- 
tropy in a fully gapped system, such as fractional quan- 
tum Hall and the gapped integer spin systemsi^ii^, in- 
dicates existence of long-range quantum entanglement 
(topological order— in equivalent parlance). Another im- 
portant progress is the demonstration that the entangle- 
ment spectrum (ES)^^^ reveals the gapless edge spectrum 
for fractional quantum Hall systems, Chern insulators, 
and topological insulatorsi^i^i2i. 

Supposing A and B to be two blocks of a large sys- 
tem in a pure quantum state, the reduced density matrix 
(RDM) PA can be obtained by tracing over degrees of 
freedom of B. Then the Von Neumann entanglement 
entropy (EE) can be computed 



Sent = -tr(p^lnpA) = -tr(/9_B Inps) 



(1) 



It has been sho'wn that for bipartite subsystems A 
and B 'with a smooth boundary, Sent has the form of 
Sent = oiL — Stop, 'where L is the length of the bound- 
ary, a is a non universal coefficient, and —Stop is a uni- 
versal constant called the topological entanglement en- 
tropy^"^'^*. Moreover, if we write the RDM in the form 
of Pa = exp(~Hent)/Z , where Z is a normalization con- 
stant, and Hent is known as the entanglement Hamilto- 
nian, the eigenvalue spectrum {si} of Hent is called the 
ES, which stores more information about the quantum 
entanglement than the EEi^. 





(a) cylinder geometry 
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(b) torus geometry 



FIG. 1: (Color online) Schematic view of a cylinder and a 
torus. The entanglement cuts divide the system into two 
equal parts A and B. For the cylinder geometry, the en- 
tanglement cut leads to one interface (a); and for the torus 
geometry, the cuts lead to two interfaces (b). 



In this paper, we study the relationship between bipar- 
tite entanglement and subsystem particle number in half- 
filled free fermion systems. It was proposed in Ref.^^, for 
systems with translational invariance in one dimension, 
the discontinuity in the subsystem particle number as a 
function of the conserved momentum indicates whether 
or not the ES has a spectral flow, which is determined 
by the topological invariant of the system^. Neverthe- 
less, this approach has an exceptional case for a half- filled 
QSH system with two-dimensional inversion symmetry. 
To overcome the inadequacy, we define spin-projected 
particle numbers, based on which spin trace indices can 
be well defined, for the QSH system with or without 
conservation. Spin trace indices are univocally related 
to the topological invariant of QSH system, i.e., the Z2 
index. We further investigate the relationship between 
the EE and subsystem particle number fluctuation. The 
latter is also dominated by the boundary excitations of 
the system, and satisfies a similar area law as the EE. 
It gives a lower-bound estimation of the EE, and can be 
measured experimentally. 

In the next section, we introduce the model Hamil- 
tonian, and explain the procedure to calculate the ES 
and EE. In Sec. HI, numerical calculation of the ES is 
carried out, and the connection between the subsystem 
spin-projected particle numbers and the topological in- 
variants in different phases is established. In Sec. IV, 
the relationship between the EE and subsystem particle 



2 



number fluctuation is discussed. The final section is a 
summary. 

II. MODEL HAMILTONIAN 

We begin with the tight-binding model Hamiltonian 
for the QSH system introduced by Kane and Mele^i^, 
plus an additional exchange field^^ 

i7 = - ^ ctcj + iVso c\(J^V^,Cj 

+ ivr ^ ct(cr X dij)zCj + ^ mictci + 5 X! ■ 

i i 

(2) 

Here, the first term is the usual nearest neighbor hop- 
ping term with ct = (Ci-^ict^) as the electron creation 
operator on site i, where the hopping integral is set to 
be unity. The second term is the intrinsic spin-orbit cou- 
pling (SOC) with v^j = {dkj X d,k)z/\{dkj x d.ik)z\, where 
k is the common nearest neighbor of i and j, and vector 
dik points from fc to i. The third term stands for the 
Rashba SOC with cr the Pauli matrix. The fourth term 
stands for a staggered sublattice potential {nii = ±m), 
and the last term represents a uniform exchange field 
with strength g. 

We consider systems with cylinder or torus boundary 
conditions, consisting of {N^ to be even) zigzag chains 
along the circumferential direction (y direction). The size 
of the sample will be denoted as TV = x Ny with Ny 
as the number of atomic sites on each chain. We per- 
form the entanglement cut along the y direction, which 
results in one or two interfaces between the two equal 
parts A and B, respectively, for the cylinder or torus 
geometry, as shown in Fig. 1. In order to examine the 
EE and ES, an Schmidt decomposition on the ground- 
state wavefunction or calculation of the RDM is usually 
needed. For non-interacting fcrmion systems, however, 
the necessary information of the entanglement can also 
be obtained from the following two-point correlators^'^ 

Cri,r2(i,j) = (cl.nCj.ra) ■ (3) 

Here, (•) means the ground-state expectation of an oper- 
ator. T can be an index of spin, pseudospin, or orbital 
degree of freedom. 

Using the Fourier transformation (FT) along the 
y direction, the Hamiltonian can be rewritten as 
H = 'Eky,i.j4ikv)'T'^,jih)'^jiky), where cl{ky) 
{cl^{ky), c] ^{ky)) are the electron creation operators. Af- 
ter performing the entanglement cut, we treat part A as 
the subsystem, and trace out the degrees of freedom of B. 
It should be noted that any of the correlators CT-i,T2(ij j) 
with i and j confined in A is unchanged by the tracing. 
When carrying out the FT on the correlators, we can get 

y h 




g 

FIG. 2: (a-c) Entanglement spectrum in the cylinder geom- 
etry (left panels) and torus geometry (right panels) for the 
QSH phase with Vso — m — 0.2, Vr = 0.1, g = (upper row), 
the insulator phase with Vr = —0.3, m — 0.3, Vao = g = 
(middle row), and the quantum anomalous Hall phase with 
Vr — g — —0.3, Vso — m — (lower row), (g) Phase diagram 
in the m versus g plane for Vso = and Vr 7^ 0. Points A and 
B correspond to the parameter values used (c,d) and (e,f), 
respectively. 

where i and j discriminate the zigzag chains. We use 
(•^I Ti (^y)''j,T"2 (^y)) to form a hermitian matrix C{ky). 
Then the entanglement Hamiltonian is given by — 

Hent = ln(C-i - 1) . (5) 

The spectrum {Q} of C is related to spectrum {ei} of 
Hent by Q = l/(e^' -I- 1), where d acts as the average 
fermion number in the entanglement energy level Si at 
"temperature" T — 1. By using the spectrum of C, the 
EE at each ky sector is given by Sent{ky) = s;, with 

= -01nO-(l-G)Ml-0) • (6) 

From the viewpoint of probability theory, Si in Eq. (6) 
can be regarded as the Shannon (information) entropy 
of the Bernoulli distribution, i.e., the i-th entanglement 
level £i has probability C^i of being occupied while (1 — Ci) 
of being unoccupied. As a result. Sent is the Shannon 
entropy of a series of such independent Bernoulli distri- 
butions. In the following, we will perform systematic 
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numerical simulations to study various phases of Hamil- 
tonian ^ in terms of the ES and the subsystem particle 
number. 



III. ENTANGLEMENT SPECTRUM AND 
SUBSYSTEM PARTICLE NUMBER 

At g — 0, Hamiltonian 1^ is the standard Kane-Mele 
model^, which is invariant under time reversal symme- 
try. The system is in a QSH phase when \m/vso\ < 
[9 — jivr/vso)"^], and is an insulator when \m/vso\ > 
[9 — j{vr/vso)'^]- On the other hand, if we set Vso = 0, Vr 
and g nonzero, a middle band gap opens when \g\ ^ \m\. 
The system is in a quantum anomalous Hall phase with 
Chern number C — ±2^^ for \g\ < |m|, and is an insula- 
tor for \g\ > \m\. The band gap closes at the transition 
point \g\ — \m\. The phase diagram for Vgo = and 
Vr ^ is plotted in Fig. [DJg). 

Figs. [2ja) and (b) show the ES for the QSH phase. 
Figs. (He) and (d) for the insulator phase, and Figs.[2fe) 
and (f) for the quantum anomalous Hall phase. Here, 
it should be emphasized that the nontrivial topological 
phases exhibit gapless ES [Figs. [2]Ja), (b), (e), and (f)], 
corresponding to physical gapless edge modes, and this 
property is named as the spectral flou^. However, the 
spectral flow is broken for the topologically trivial phase 
[Figs. [SJc) and (d)], which is also consistent with the 
property of the correspondent edge states. 

In a recent work^-, the authors proposed a new char- 
acteristic quantity called the "trace index" to describe 
topological invariants, which is defined through a subsys- 
tem particle number operator NA{ky) = I^igA 4,fc„Ci,fcj, . 
The expectation of NA{ky) is given by 

(NAiky)) = {GS\J2cUky)c,iky)\GS) - TrC . (7) 

In Fig. [31 we plot the expectation of NA{ky) for the 
three different phases mentioned above. In the cylinder 
geometry, NA{ky) is discontinuous at some discrete mo- 
menta in the nontrivial topological phases, as shown in 
Figs.lSfa) and (c). This is in contrast to the normal insu- 
lator phase [see Fig. [3ljb)], where NA{ky) is a continuous 
function of ky. In the torus geometry, NA{ky) is exactly 
equal to half of the total particle number in the ky sector, 
without showing any discontinuity, because the change of 
the particle number in A around interface / is just can- 
celed by that around interface II due to the rotation in- 
variance of the torus. In the cylinder geometry, the trace 
index was defined as the total discontinuities of {NA{ky)) 
with varying momentum. Alexandradinata, Hughes, and 
Bernevig^^ presented a detailed analysis and proved that 
the trace index is equivalent to the Chern number (or 
Z2 invariant) for the Chern {Z2) insulators. Therefore, 
the subsystem particle number provides a new alternative 
tool to reveal the topological invariants. 

However, as mentioned in Ref.— , there is an excep- 
tional case in which the subspace of the occupied bands 




FIG. 3: (Color online) The subsystem particle number in the 
cylinder geometry and torus geometry for the QSH phase (a), 
the insulator phase (b), and the quantum anomalous Hall 
phase (c) with all the parameters same as those in Fig. 2. (d) 
is also for the QSH phase with Vso = Vr — 0.2, where the two- 
dimensional inversion symmetry is retained. Discontinuities 
in the expectation of the particle number can be observed 
only in the cylinder geometry. 



at the symmetric momenta is not closed under time re- 
versal in the ground state. If at the symmetric momenta 
the Kramers' doublet that extends along the edge of A 
is singly-occupied, {NA{ky)) is continuous, even when 
the system is in a nontrivial topological phase. For the 
half-filled system under consideration, an exception still 
happens. While the two-dimensional inversion symmetry 
remains unchanged (to = 0), NA{ky) becomes continu- 
ous, as shown in Fig. |31[d). This is because the Kramers' 
partners extending along the edge simultaneously cross 
the Fermi level at the symmetric momentum {ky = tt) 
and have opposite contributions to the discontinuities of 

{NA{ky)). 

To overcome this difficulty, enlightened by the defini- 
tion of the spin Chern number—, we define spin trace 
indices. Considering that Sz is not necessarily conserved, 
it is an adaptable way that we choose operator PszP 
to split the fiber bundle of the occupied states into two 
bundles with nontrivial Chern numbers, where P is the 
ground state projector. At half filling and in the pres- 
ence of time reversal symmetry (3 = 0), PszP is always 
a time-odd operator {TPszPT~^ — —PszP), so that the 
spectrum of PszP is symmetric in respect to the origin. 
As a result, the spectral spaces of PszP can provide a 
splitting of the Hilbert space spanned by the wave func- 
tions of the occupied states, resulting in a smooth de- 
composition P{ky) into 



P{ky)=P+{ky)®P-{ky) 



(8) 



with a — ± corresponding to the positive and negative 
sectors of the spectrum of PszP for all ky S (0, 27r]. 
Straightforwardly, the two-point correlator matrix can 
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FIG. 4: (Color online) Subsystem spin-projected particle 
numbers in the cylinder geometry for the QSH phase (a- 
c) with g = and different parameters: (a) Vso = 0.2 
and Vr = m — 0, (b) Vao = Vr = 0.2 and m = 0, (c) 
Vso — m — 0.2 and Vr — 0.1, and for the insulator phase 
(d) with Vso — Vr = 0.05, m = 0.5 and g = 0. 



also be decomposed into C{ky) — C^{ky) ®C^{ky). In 
this way, the Chern numbers C± of the two spin sec- 
tors can be well defined on these new wave functions. It 
has been proved that C± are topological invariants pro- 
tected by the energy and spin spectrum gap&ii. It will 
be shown below that the traces of C^, called the spin- 
projected subsystem particle numbers, are related to the 
topological invariants. We plot TrC" (a = ±) as func- 
tions of ky in Fig. 01 Both TrC"(fcj,) are discontinuous in 
the QSH phase [Figs. lUJa-c)], in contrast to the contin- 
uous functions in the normal insulator phase [Fig. Ul^d)]. 
If TrC"(fcj^) is discontinuous at some discrete momenta 
{kdis} with kdis G (0, 27r], we can define the spin trace 
indices as the total discontinuity, i.e., difference between 
the limits of TrC"(fcy) from right and left. 




in the thermodynamic limit. As shown in Figs.lH^a) and 
(b), no matter whether Sz is conserved, both TrC+(fcj,) 
and TrC~(fcy) show discontinuities at momentum ky = t: 
with = 1 and A~ = — 1 in the QSH phase, where the 
two-dimensional inversion symmetry is present (m = 0). 
Figure HJc) shows the discontinuities of TrC"'"(fcj,) and 
TtC~ (ky) in the QSH phase in which Sz is not conserved 
(vr ^ 0) and the two-dimensional inversion symmetry is 
broken (to 7^ 0). In this case, the spin trace indices are 
equal to 1 and —1, respectively, contributed by two differ- 
ent momentum points. It is noteworthy that in analogy 
with the Laughlin gauge experiment. A" can be regarded 
as the number of particles with spin a pumped from one 
edge to the other when a unit fiux is inserted adiabati- 
cally, and so A"' is equivalent to the Chern numbers Cq. 



On the other hand, the Z2 index can be defined as the 
parity of A", 

Az^ = A°'mod 2, (10) 

for any a. As shown in Figs. SI Az2 — 1 for QSH phase 
[Figs.m^a-c)] and Az^ = for insulator phase [Figs.|4l[d)]. 
Therefore, the subsystem particle number expectation 
can be used to characterize the topological invariants. 
Especially, for the QSH systems, the spin trace indices 
are well-defined quantities that can reveal the Z2 invari- 
ant and distinguish different quantum phases. 

IV. ENTANGLEMENT ENTROPY AND 
SUBSYSTEM PARTICLE NUMBER 
FLUCTUATION 

We have shown that topological properties of the 
ground state can be extracted from the expectation of 
subsystem particle number. Now we turn to the variance 
of NA{ky). In the past two years, extensive works have 
been devoted to the study of the relation between the EE 
and subsystem particle fiuctuation for non-topological 
systemE^'5. In this section, we will show that the rela- 
tion is rather general, it does apply to non-interacting 
electron systems with a nontrivial band topology. We 
start from the definition of the variance 

AiVKfc,) = {Nl{ky)) - {NA{ky)f . (11) 

Substituting Eq. ^ into Eq. dTT]) and using the Wick's 
theorem to expand all the four-point correlators, one can 
obtain 

= Tr[C(l-C)], (12) 

yielding AN\{ky) = J^idi^ - 0), which is in keeping 
with the variance formula of the Bernoulli distributions. 
It then follows that the variance is also dominated by the 
low-energy boundary excitations (0 < Ci < !)• Moreover, 
each maximally entangled state with e„i — (Cm = 1/2) 
contributes a maximal value to the subsystem particle 
number fluctuation and the EE, which cannot be elimi- 
nated by adiabatic continuous deformation. 

In order to find a definite relationship between the EE 
and the variance, one can construct a concave function 
f{x) — — lna;/(l— x) for x e [0, 1], and apply the Jensen's 
inequality 

- a;lna; - (1 - a:)ln(l - x) ^ (41n2) • x{l - x) . (13) 

The equality is taken if and only if x — 1/2. Equation 
enables us to make a lower-bound estimation of the 

EE 

Sentiky) 1^ {iln2) ■ ANliky) . (14) 
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FIG. 5: (Color online) Entanglement entropy in compari- 
son with subsystem particle number fluctuation for the QSH 
phase (upper row) , the insulator phase (middle row) , and the 
quantum anomalous Hall phase (lower row), in the cylinder 
geometry (left panels) and torus geometry (right panels). All 
the parameters are same as those in Fig. 2. 



ial system. 

Furthermore, one can use NA{ky) = J2i(=A'4k 
verify AN\ = Efc, A7Vl(fc,) ^ ^ / dfc,7Vl(fc^), indi- 
eating that /\N\{ky) satisfies a area law2i, similar to the 
EE, Sent = Y^ky Sent{ky) ^ J dkySent{ky). Therefore, 
the subsystem particle number fluctuation shares several 
common characteristics with the EE, and so can be uti- 
lized to detect the EE experimentally. 



V. SUMMARY 

To conclude, we have investigated the relationship be- 
tween the quantum entanglement and subsystem particle 
number. The spin trace indices can reveal the topological 
invariants and be used to classify different phases in QSH 
systems. This new tool always works well even though Sz 
is not conserved. As to the subsystem particle number 
fluctuation, it shares several common properties with the 
EE. They both satisfy the same area law, and are domi- 
nated by the boundary excitations with each zero mode 
having a maximal contribution. The connection between 
the two quantities is universal, regardless of whether the 
system has a nontrivial band topology. As a result, the 
subsystem particle number fluctuation, as an observable 
quantity, can be used to detect the EE experimentally^^. 



Thus a lower bound of the EE is given by sa{ky) = 
(4 In 2) • AN'^{ky), which is directly proportional to the 
particle number fluctuation of subsystem. In Fig. [5] we 
plot Sentiky) and so{ky) in the QSH phase, insulator 
phase, and quantum anomalous Hall phase. In all the 
cases, the curves for the particle number fluctuation be- 
have somewhat similarly, and are very close to the cor- 
responding EE. This similarity was observed in the non- 
topological systems lately^^, and here we find that the 
similarity remains to hold for the topologically nontriv- 
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